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Abstract 

> 

' The branching ratio is calculated for three different models of 2d gravity, using 

dynamical planar (p graphs. These models are pure gravity, the D = —2 Gaussian 
■ model coupled to gravity and the single spin Ising model coupled to gravity. The 

, ratio gives a measure of how branched the graphs dominating the partition function 

I are. Hence it can be used to estimate the location of the branched polymer phase for 

the multiple Ising model coupled to 2d gravity. 

Oh! 

^ ■ 1 Introduction 



X 



Many interesting results in the area of 2d quantum gravity have been obtained by studying 



H ' models of matter fields coupled to dynamical triangulations. However, all the exactly 

solved models have had central charges, C, of less than or equal to one (if one ignores 
models for which C = oo [|^, |2|). Using conformal field theory it is possible to make 
predictions for the critical exponents in models of matter coupled to 2d gravity ||3|, ^, 
the so-called "KPZ formulae", but again these formulae are only valid for C < 1. The 
nature of this "C = 1 barrier" remains something of a mystery. Many models of matter 
coupled to gravity have a tendency to degenerate into branched polymer configurations 
for large enough C. It is thus tempting to hypothesise that a branched polymer phase 
begins at C = 1, causing the break down of the KPZ formulae here. Models in which 
multiple Ising spins are coupled to dynamical triangulations (or their dual cf)^ graphs) 
provide a convenient way of testing this hypothesis; such a model with p independent 
Ising spins on each vertex of the graph, has a central charge of C = p/2 associated with 
it. However, Monte Carlo simulations for values of p up to p = 16 have failed to show 
any convincing evidence of the existence of such a branched polymer ("BP") phase. Since 
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there are reasons to believe that such a phase exists [||, there are two possibihties to 
consider: firstly, that there exists an intermediate region between C = 1 and the start of 
the branched polymer phase, or secondly, that the BP phase does indeed start at C = 1 
and that the MC simulations have for some reason failed to detect it. This paper aims to 
throw some light on the question by examining the branchedness of the graphs for various 
exactly solvable models. 

In a recent paper we studied correlation functions in a model of multiple Ising 
spins on a dynamical planar c/)^ graph. In the paper, a function B was defined which 
played a key role in the calculation of the critical exponents; B gives a measure of how 
branched the 0^ graphs contributing to the partition function are and hereafter it will be 
called the "branching ratio". The branching ratio lies in the range to 1, and throughout 
the branched polymer phase B equals one, when evaluated at the critical value of the 
cosmological constant. Thus the branching ratio gives a measure of how close to the BP 
phase we are, and by evaluating it for various models we will gain some indication of the 
extent of the branched polymer phase for the multiple Ising model coupled to gravity. 

In this paper we calculate the branching ratio, which is defined in the next section, for 
three different models (all of which are defined in terms of planar graphs): pure gravity 
(section ^), the D = —2 Gaussian model coupled to gravity (section ^ and the single spin 
Ising model coupled to gravity (section |5|). As might be expected the branching ratio 
increases with increasing central charge. In section ^ we conclude by using these results 
to make a prediction of the location of the branched polymer phase for the multiple Ising 
model coupled to gravity; it is estimated that there exists a BP phase for central charges 
above some value C, where C lies in the range 1 < C <,2.75. 

2 The branching ratio 

In our previous paper Q we used a simplified definition of distance rather than the usual 
geodesic distance. The geodesic distance between two vertices A and B oi a cj)^ graph is 
defined to be the shortest path along links, from ^ to -B, counting one unit of distance 
for each link traversed. However, in our definition of distance, only links separating two 
one-particle irreducible ("IPI") subgraphs are counted; figure ^ gives an example, each 
shaded circle represents a IPI subgraph, hereafter referred to as a "blob". 

Using this definition of distance it was possible to show that the number of blobs at 
a distance r from the root, denoted (n)^, was given by (n)^ = B^' , where i? is a function 
of the various parameters in the model. 

For the case of a single Ising spin coupled to each vertex of a graph (using the set, Q^, 
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Figure 1: Measuring distances from the root blob 
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of all rooted planar (fi graphs), we derived the following formula for B 
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1 - 3xr + 
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Jfi' 



(1) 



where T = Z^{x,/3). The grand canonical partition function for the single spin Ising 
model on the set of graphs is defined by 



Gagi 



(2) 



{5} <ij> 



where Si is the spin on vertex i and t = tanh This will be referred to as model I. The 
coupling constant /3 is the inverse temperature of the model and each vertex in the graph 
is weighted with a factor of x. The first summation is over graphs G in the set Q^, where 
N is the number of vertices in G, the second summation is over all spin configurations 
and the product is over nearest neighbour pairs. Note that there is no spin at the end of 
the root and hence the root link is not included in this product. The partition function 
Z^^ = Z^^{x',f]') is defined as in (|2|), but using the set of planar rooted IPI cf)^ graphs, 
Qr', this will be called model II. The renormalized coupling constants t' = tanh/?' and x' 
are given by 
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x{l - 2xT)-2. 



(1 - 2xTt) 

and the following equation which relates Z^ and Z^^ was derived in our paper Q, 

T = Vl - 2xT Z"{x', 13') + xT'^. 



(3) 



(4) 
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3 Pure gravity 



In the case of pure gravity (that is, equation (^) with /3 = 0) equation ([T|) reduces to 



B{x) 



xT 



1 - 2xT 



1 - 3xT + 3(1 - xT) 



dZ 



II 



Z" dx' 



(5) 



We are interested in evaluating this at the critical value of x for model I, which is denoted 
by xi- Prom reference [|, = 1/(2.33) and T{xi) = 33 (l - ^3/2) . The corresponding 

value of x' is x^^ = and in reference it was shown that 



x' dZ" 



5 
3' 



(6) 



Z^^ dx' 

Thus writing Be for B{x],), we have that the branching ratio at the critical value of x 
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5c = 1 - ^ ~ 0.4226 
V3 



(7) 



for pure gravity, with rooted planar i;^'^ graphs; this model, of course, corresponds to a 
central charge of zero. 



4 Gaussian model 



In this section we will derive B^. for the Gaussian model in Z? = —2 dimensions, which 
has a central charge of —2. For this choice of D the partition function has an especially 
simple form Q, 

Again the set of all rooted planar c/)^ graphs is being used and in addition there is a 
summation over the trees T in the set of all trees spanning G; this set is denoted by Tq- 

The partition function has an integral formulation |1C], 



1 /"^ 

Z^(x) = — / dtVl-t"^ f 1 - Vl - 8xt) , 

^ XTT 7_1 V / 

which has a critical value of x given by = |. At this critical value. 
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Figure 2: T = f{T) 
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and diverges. 

A formula for B{x) in this model can be obtained by using a derivation similar to 
that for the multiple Ising model Define a function, f{y), which takes a rooted IPI 
blob (both graphs and spanning trees are being summed over in the blob) and glues an 
arbitrary number of legs, each weighted with a factor y, on to the blob. Again we define 
T = Zg{x) to save writing and have that T = f{T), which is illustrated in figure |2|. To 
derive a formula for /(y) we need to consider the effect of adding an arbitrary number of 
legs to a link. For a given graph G and spanning tree T there are two different types of 
link, those that form part of the spanning tree, which will be referred to as "black" links 
(drawn as thick lines), and those which are not part of the spanning tree (these are called 
"white" and drawn as thin lines). Adding legs to a black link changes the contribution it 
makes to the partition function (fig 0a) : 



1 ^ 1 + 2xy + {2xyf + 



1 



1 — 2xy 



(12) 



For each leg added there is a weight of y, a weight of x for the extra vertex and a factor 
of 2, since the leg can be hung from the link in one of two directions; this accounts for 
the factors of 2xy. 

Similarly for white links there is the change (fig Hb): 



1^1 + 2{2xy) + 3{2xyf + 4(2xy)^ H 



^1 - 2xy)2' 



(13) 



In this case one must also take into account the different ways of connecting the new legs 
to the existing spanning tree. 

The partition function for rooted IPX blobs is defined by 



.N 



(14) 



Adding legs to this changes the partition function for the blob to 

J2 E x^(l-2xy)-(^^+2^^), 
Gegii TeTc 



(15) 
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Figure 3: Renormalization of links: (a) Black, (b) White 



(a) 



















+ - 




- 


+ _ 


+ 


_ + 

















+ 



(b) 

















> 


+ - 


+ - - 


f - + 






+ _ 



where the number of black links is = — 1 and the number of white links is Lw 
i(iV + 1). Thus Lb + 2Lw = 2N and hence 



f{y)=xy^ + Z'J (^x{l-2xyy 



(16) 



It is convenient to define x' = x(l — 2xT) ^ . Note that at x = x^, we have x' = x^^ 

2 



, which agrees with the result in reference |11]. As before B is given by 



dy 



y=T 



(17) 



SO that, with = Z^^x'), 



B{x) = 2x 



' dZ" 2x' 

T + — ^- 



dx' l-2xT 



(18) 



Differentiating the equation T = f{T) with respect to x' yields an expression for -jjp- and 

0.1670, (19) 



then evaluating (18) at x = gives 



Bc = l 



3n_ 
8V2 



for the D = —2 Gaussian model. 



5 Ising model 

Finally, we will calculate the branching ratio for the Ising model on planar cp^ graphs 
with a single spin per vertex; the matter in this model has a central charge of ^. To use 
equation (|^) we need the partition functions for both model I and model II, which we 
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obtain by solving the following two-matrix model, 



VA VB exp TV 



--{A^ + B^) + cAB 



+ l-^lA^ + B'^^+xVMiA + B) 
3^ V ) 



(20) 



where A and B are M x AI Hermitian matrices and we wish to take the planar limit in 
which M — > cxD. Each vertex is weighted with g and the coupling constant c equals e~^^. 
In the case X = 0, F will give the model I partition function up to various factors; this 
has been solved in references 13 1. Alternatively by choosing A so that graphs with 
tadpoles are cancelled, we gain the partition function for model II. 

This matrix model can be solved using orthogonal polynomials [Q. In the planar 
limit and dropping the constant term, we have 

.1 n 



Fig,c,X) = hn 



1 



9 

1/1^2 



-u^ — u 



— —du H h - In 1 - c"' 



3 \2g 



where 



ug 
1 

8p 



(1 + po) [{I - c) (1 + po) - 8Xg 
2 - z(/>2 - c^) + 2cz^ 



(p + 1) (p - 1 + 2c) + 4Xg 



and Po is p evaluated at u = (i.e. at z = 0), 

Po = -c- ^{l-cf -AXg. 
After some considerable algebra we obtain. 
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dF\ 
dx), 



1 + Pi - -^zfpiiPi 



(21) 

(22) 
(23) 

(24) 
(25) 



g,c 9^ 

where pi = p{u = 1) and z\ = z{u = 1). Note that F generates unrooted (i.e. vacuum) 
graphs and that by differentiating with respect to A we generate rooted graphs instead. 

Now, the model I partition function that we want to calculate is related to that 
generated from the matrix model by 



9 



OF 
9A 



A=0 



Thus dropping the subscripts, 

1(1 + /') 



xT 



2(l-c) 



Pip-c)ip-a) 
(4/9^ + c/j^ — 2c^p — ca) 



(26) 



(27) 
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where c = e a = 1 — 2c, 



z= ^(/9 + l)(p-a), 



/ = z{p^ - (?) + 2cz^ = 2(1 - cfx^. 
At x^, which is given by =0, 



^3 , o„„2 



3(1-20 



2/>^ + 3cp^ + c(l - 2c) = 
There is a third order phase transition in model I, at 



(for /3 > /?*), 
(for /3 < /?*). 



c* = — f 2^7 - 1 
27 V 



between a magnetized and an unmagnetized phase. 
For model II we need to choose A = A such that ^ 



vanishes. This gives us, 



p(l+p)(p + c)2(p-c) = 



and 



mA- + 



4 (p2 _ 2pc - c) 
At the critical value of x, denoted x^^, 

+ 2p + c = 
- 3c/9^ - 3cp + c^ = 



2(p2 _ 2/jc - c)2 

3/9^ + - pc- c(l - 2c)) . 



(for /? > /?*), 
(for /? < 



There is a third order phase transition at c* = and the model II partition funct 



(l-c) 
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Defining 



E{x,P) 



1 



3x 



(dz"\ 



+ 2.(t-l)(^^| 



^^^(x,^) 

we have from equations (|3^), (|^) and (^7|) that 

4(l-c)(3p3 + c2) 



E{x,(3) = 4c - 3 + 



(3p3 + p2 _ pc- c(l - 2c))' 



Then 



5(x,/3) 



xT 

1 - 2xr 



1 - 3xr + (1 - xT)E{x',l3' 



where E is given by the previous equation, x' and /?' by @ and xT is given by (27). 
Evaluating this at the critical value of x gives Bc{(3) = B{x^^{(3), (3); figure ^ gives a plot 
of Be against c = e~'^^. The branching ratio is a maximum at the phase transition and 
one can show that the branching ratio there is 



Bcif3* 



1 



4/7 



0.5275. 



Note that as expected both model I and model II have ^jstr 
and 7*^r = — | at the critical point. 



(41) 



i away from the transition 



Figure 4: Graph of the branching ratio, Bc{P), against c = e for the single spin Ising 
model coupled to 2d gravity. 
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6 Conclusion 



In this paper we have calculated the branching ratio for three different models coupled to 
2d gravity, in the form of rooted planar cj)^ graphs. In figure |5| are plotted the values of B^. 
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Figure 5: Graph of the branching ratio, B,,, against the central charge, C, for various 
models coupled to 2d gravity. 
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against the central charge, C; for the case of the single spin Ising model we used Bc{(3*), 
since the graphs are most branched at the phase transition. We see that the graphs 
tend to become more branched as the central charge increases. This is not unexpected 
since many models coupled to gravity degenerate to branched polymers for large enough 
central charge. One should perhaps note at this point that Be is not a universal function 
of the central charge, so that different models with C = —2 will in general have different 
branching ratios. Thus the data point for the Gaussian model is probably not very reliable 
when trying to determine the behaviour of the multiple Ising model. For this reason we 
extrapolate using only the last two data points and estimate that Be = 1 at C = 2.75 (the 
extrapolation is the right-most dotted line in figure ||). This is likely to be an overestimate 
since the graph seems to be curving upwards with C. Even if we ignore the data point for 
the Gaussian model we know that for C — > —oo, Be must tend to a finite non-negative 
value and in this limit one might guess that Be 0; so that it seems likely that the graph 
will be curving upwards in the region of interest. In our previous paper Q it was shown 
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that regions of the phase diagram with Be = 1 correspond to the branched polymer phase 
(7str = f ) and its boundary {'jstr > 0). It is known that there is no branched polymer 
phase for C < 1 and thus the left-most dotted line indicates an approximate upper limit 
on possible values of B^- So we conclude that the curve of B^ against C probably lies 
within the region bounded by the two dotted lines and hence that there exists a branched 
polymer phase for central charges greater than some value C where 1 < C ^2.75. Note 
that this estimate for C is lower than previous cruder estimates 15|. Our result is 



still consistent with the hypothesis C = 1, although it is not possible to rule out the 
existence of an intermediate region. Even though we are unable to distinguish between 
these two possibilities it seems likely that any intermediate region which exists would be 
quite narrow. Further work will be required in order to determine more fully the phase 
diagram for the multiple Ising model coupled to gravity and to understand the nature of 
the C = 1 barrier in general. 
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